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Abstract. We present a large-scale study of the IC 1396 region using new deep NIR and optical images, complemented by 
2MASS data. For ten globules in IC 1396 we determine (H-K, J-H) colour-colour diagrams and identify the young stellar 
population. Five of these globules contain a rich population of reddened objects, most of them probably young stellar objects. 
Two new HH objects (HH 865 and HH 864) could be identified by means of [SII] emission, one of them a parsec-scale flow. 
Using star counts based on 2MASS data we create an extinction map of the whole region. This map is used to identify 25 
globules and to estimate their mass. The globule masses show a significant increase with the distance from the exciting 06.5V 
star HD 206267. We explain this correlation by the enhanced radiation pressure close to this star, leading to evaporation of the 
nearby clouds and hence smaller globule masses. We see evidence that the radiation from HD 206267 has a major impact on 
the star formation activity in these globules. 
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1. Introduction 



formation - stars: winds, outflows - ISM: Herbig-Haro objects - ISM: jets and outflows — Individual: 



I1979> . Large-scale observations of this region in the rotational 



Star formation takes place not only in giant molecular clouds 
but also in small isolated globules. Works e.g. by Sugitani et 
Schwartz et al 



al. J199U . Schwartz et al. J199U and others showed that such 
places are associated with young stellar objects. The identifica- 
tion of young stellar clusters or outflow activity in such glob- 
ules hints an ongoing star formation process. Star formation in 
globules might be induced by the propagation of an ionisation 
shock front, the so-called radiation driven implosion mecha- 
nism (Reipurth ll983> . However, it is not fully clear what the 
main properties are that influence star formation within these 
globules (e.g. density, mass, or size of the globule, strength of 
the ionisation shock front). Further, it would be interesting to 
see if and how these properties influence the number and clus- 
tering of the forming stars. The investigation of a larger, homo- 
geneous, and as far as possible unbiased sample of globules is 
an ideal way to obtain a deeper understanding of these issues. 

IC 1396 is one of the youngest and most active HII re- 
gions in the Cep OB 2 group of loosely clustered OB stars 
(Schwartz et al. 11991b . The nebula is excited by the 06.5V 
star HD 206267 (Walborn & Panek [19841 and contains 15 
small clouds and globules associated with red IRAS sources 



(Schwartz et al. U99U at a distance of about 750 pc (Matthews 
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CO lines were done by Patel et al. (19951 and Weikard et al. 
( 1996). The numerous sharp-rimmed clouds and the relative 
proximity make this region an ideal place to study star forma- 
tion in a large, homogeneous sample of small globules. 

Some of the globules in IC 1396 were already investigated 
in detail (IC 1 396N: Beltran et al. I2U021 Nisini et al. |2"0TJT1 
Codella et al. 120001 IC1396W: Froebrich & Scholz 120031 
IC 1396 A or the Elephant Trunk Nebula: Nakano et al. IT9891 
Hessman et al. 119951 Reach et al. 2004 1 and/or are known 
to harbour outflow sources (IRAS 21388+5622, Duvert et al. 
IT990I Sugitani et al. IT9971 De Vries et al. 120021 O gura et al. 
120021 IRAS 22051+5848, Reipurth & Ballv l200!V Relatively 
little is known about most of the other globules, which are 
therefore targets for this new study. As shown e.g. in Froebrich 
& Scholz (2003), deep NIR imaging in JHK and the construc- 
tion of (H-K, J-H) colour-colour diagrams can reveal embedded 
young objects or clusters in such globules. These data can be 
used in conjunction with star counts (see e.g. Kiss et al. [2000 1 
to estimate the extinction and hence determine the mass of the 
globules. Simultaneous observations of the region in narrow- 
band filters centred either on the 1-0 S(l) line of molecular hy- 
drogen at 2.122 fim or the [SII] lines at 671.6 and 673.1 nm will 
uncover outflow activity from young stars. 

This paper is structured as follows: In Sect.|5]we describe 
our data obtained in the optical and NIR and specify the data 
reduction process. Our photometry in the NIR images with 
emphasis on (H-K, J-H) colour-colour diagrams is shown in 
Sect. [3] The creation of NIR extinction maps to estimate the 
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mass of the globules is described in Sect.0] followed by a dis- 
cussion of the determined globule properties in Sect.[5]and the 
characterisation of the newly discovered outflows in Sect.[6] 
Finally, the results are discussed and summarised in Sect.Q 
Some more detailed technical descriptions of the data analysis 
procedures can be found in AppendiceslAl-ICl 

2. Observations and data reduction 

2. 1. Near infrared data 

The near infrared (NIR) observations were made towards 10 
out of the 15 IRAS sources listed in Schwartz et al. (I199U . 
Two of them have been studied before (IC 1396N, Beltran et 
al. 120021 Nisini et al. I200T1 Codella et al. 120001 IC 1396 W 
Froebrich & Scholz 120031 . and another three globules could 
not be observed due to bad weather conditions during our run. 
The observed globules are listed on the top of Tabled One 
main purpose of the NIR data was to construct colour-colour 
diagrams for each globule in which we want to detect reddened 
sources. This can be done by comparing the measured colours 
of the targets with that of the main sequence. To get an esti- 
mate for the colours of the main sequence, we observed main 
sequence stars with known spectral type in the same way as our 
IC 1396 targets. These standard stars were selected from the 
SIMBAD database and deliver an unreddened main sequence 
for our colour-colour diagrams. We preferred to define this 
main sequence by observations rather than with theoretical evo- 
lutionary tracks, because this avoids mismatches from inconsis- 
tent photometric systems, which can be substantial, as we have 
shown in Froebrich & Scholz d2003> . We selected the standard 
stars so that they are generally less than 10 degrees away from 
the IC 1396 region. Thus, they can be observed roughly at the 
same airmass as IC 1396, which excludes systematic offsets be- 
cause of differential extinction. 

Our near infrared data were obtained on six nights from the 
18th to 23rd of July in 2003 with the 2.2-m telescope on Calar 
Alto, Spain. We observed with the MAGIC camera (Herbst et 
al. [19911 in its wide-field mode (6.'92 x 6.'92 FoV). Using a 
3x3 dither pattern around the central coordinates with a shift 
of half a detector size we obtained 13.'5 x 13.'5 sized mosaics 
of the IC 1396 globules. This field size is comparable with the 
typical diameter of small globules in this region (see Tabled- 

For the standard star observations, we chose a smaller shift 
between the images to reduce overhead times. All fields were 
observed in J, H, and K', the IC 1396 globules additionally in 
a narrow -band filter centred on the 1-0 S(l) line of molecular 
hydrogen at 2.122 fim (hereafter called the F^-filter). The K' 
filter is very similar to the 2MASS Ks band and will be called 
K in the following. The total per-pixel integration time in the 
IC 1396 mosaics was 324 s in each broad-band filter. The inte- 
gration times for the H2 -filter were at least 2160 s per pixel, in 
most cases we reached 3240 s per pixel. For the standard stars, 
appropriate integration times were chosen to avoid detector sat- 
uration. 

The weather conditions during the observations of the stan- 
dard stars were photometric. Our ten globules in IC 1396 were 
also mainly observed under photometric conditions. For each 



globule, we obtained at least one full set of mosaics (JHK), and 
we took care to assure that at least one of these sets was ob- 
served under photometric conditions, so that a calibration of the 
whole dataset is possible. In case of cirrus, some of the globules 
were observed longer to obtain a uniform limiting magnitude in 
all fields. The narrow-band observations of the globules in the 
H2 -filter were done mostly under non-photometric conditions 
or at higher air mass. The seeing conditions during our run were 
excellent (below 1 "), but we are limited by the large pixel scale 
of 1'.'6 per pixel. All broad-band images, IC 1396 globules as 
well as standard stars, were obtained in the airmass range be- 
tween 1.07 and 1.5. The positions of the observed globules are 
shown in Fig. ^( small squares). 

Our standard NIR data reduction included flat-fielding (us- 
ing sky-flats), sky-subtraction, and co-addition to mosaics. To 
co-center the images into a mosaic we used all stars in the 
field to ensure a high astrometric accuracy. Mosaicing and sky- 
subtraction was done using the xdimsum package in IRAF 1 
(Stanford et al. I1995i . For each filter and globule we produced 
mosaics of the 3x3 dither patterns. Additionally, all single mo- 
saics of each filter were stacked to a deep mosaic. A plate so- 
lution for each mosaic was obtained using the 2MASS sources 
in the field. 

2.2. Optical data 

A large-scale optical survey of the IC1396 region was car- 
ried out with the Schmidt camera at the 2-m telescope of the 
Thiiringer Landessternwarte Tautenburg in two observing runs 
(9.-12. Sept. 1999 and 21.-25. June 2004). We observed with 
a narrow-band filter centred on the [SII] emission line as well 
as with a broad-band I-filter, to be able to identify emission 
regions. These optical observations cover about 12sqdeg. The 
survey field is shown in Fig.^ A standard data reduction was 
performed including bias subtraction and flat-field correction. 
Emission features were searched for by comparing the narrow- 
band [SII] images with the I-band images. 

3. Photometry 

3.1. Technique 

For source detection, we applied the SExtractor (Bertin & 
Arnouts 1996 1 to our deep K-band mosaic of the globule. The 
relative offsets between the positions in this 'master' mosaic 
and all other mosaics of the same field were determined by 
measuring the pixel coordinates of a bright star in each image. 
Applying these offsets to the source catalogue of the 'master' 
mosaic, we obtained catalogues for all mosaics of the field. 
Subsequently, we performed aperture photometry for each ob- 
ject in the catalogues using the daophot package within IRAF. 
Because of the large pixel scale, the seeing is more or less con- 
stant in all mosaics. Therefore, we decided to hold the aperture 
for the photometry constant for all mosaics. Sky coordinates 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 1. IC 1396 field investigated in this paper. The background image is the extinction map obtained by accumulated star counts 
using the 2MASS J-band data (see Sect.|4j. Black represents regions with high extinction, gray zero extinction, and white negative 
extinction (caused by star clusters). The small squares mark the fields observed with MAGIC in JHK and in the 1-0 S(l) line of 
H2 (IC 1396 W, the western most field, is marked also, since it is included in the discussion). Circles correspond to the globules 
listed in the bottom half of Table^(a few globules are slightly outside the region shown here). The numbers at the side of circles 
and small squares are identical to the globule numbers used in Tabled The area marked with dashed lines was observed with 
the Schmidt telescope in Tautenburg in I and [SII] . The very dark, circular region at about (2l'?73:58.8°) is not a globule but the 
very bright star /1 Cep, which caused non-detections of stars in its vicinity. A similar fake globule is generated by the star £ Cep 
at (22 1 ?19:58.2°). The + signs mark all IRAS sources in the field. 



for the sources were determined by applying the plate solu- 
tion to the pixel coordinates. As the result we obtained for each 
globule a database with sky coordinates and instrumental mag- 
nitudes in JHK and H2. 

The instrumental magnitudes of the standard stars were 
measured in a very similar way as the IC 1396 object values, 
i.e. with aperture photometry using daophot. Since all these 
stars were observed under photometric conditions and airmass 
similar to the IC 1396 object fields, their colours can be directly 
compared with them. According to the spectral types given in 
the literature, our standard stars span a spectral range from BO 
to K8, hence they can be used to establish a main sequence in 
the colour-colour diagram. 

Tests were performed to verify the reliability of our instru- 
mental photometry (e.g. repeated observations, field overlap). 
Only in one of the ten fields (IRAS 21539+5821) were larger 
systematic photometry mismatches found, and hence the field 
was excluded from further analysis. For a more detailed discus- 
sion see AppendixlAl 

3.2. Colour-colour diagrams and reddened sources 

We show colour-colour diagrams for five of the nine globules, 
which contain more than 20 reddened objects, in Figs. 12161 
These figures contain a) datapoints for the IC 1396 objects as 



small dots, b) the main sequence datapoints as filled triangles, 
c) the extinction path calculated with the extinction law given 
by Mathis 1 1990) as dotted lines, d) reddened objects as filled 
squares and arrows. To allow a reliable comparison between the 
objects in IC 1396 and the standard stars, these diagrams are 
shown in instrumental magnitudes, which avoids colour mis- 
matches. A rough transformation into the 2MASS photomet- 
ric system can be obtained by adding about 0.3 mag to the H- 
K values. The diagrams show all datapoints in the 13.'5x 13.'5 
field for each globule. We generally selected objects with H- 
K > 0.5 mag as reddened sources. A more detailed discussion 
of the colour-colour diagrams and transformation of instrumen- 
tal to 2MASS colours can be found in AppendixlBl 

As noted above, we found a large population of reddened 
objects for five of the nine globules. The remaining four glob- 
ules contain only a small number of reddened sources. To as- 
sess the nature of all these red objects, we estimated the number 
of field stars in the direction of IC 1396 using the Besancon 
Galaxy model of Robin et al. d2003l >. This model is avail- 
able online 2 and is able to generate a photometry catalogue 
of objects in a certain Galactic direction with user-defined 
colour constraints. Our selection criterion for red objects was 
H-K > 0.5 mag in instrumental colours, which translates into 
about H-K > 0.8 mag in absolute colours. 

2 http://www.obs-besancon.fr/www/modele/modele.ang.html 
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Fig. 2. Colour-colour diagram of the IRAS 2 1 324+57 1 6 glob- 
ule. Small dots are objects in and around the globule. Filled tri- 
angles are field stars with known spectral type (which mark the 
unreddened main sequence). Filled squares, arrows up (lower 
limit in J-H), and arrows right (lower limit in H-K and arbi- 
trary J-H) mark reddened objects. The dashed line shows the 
criterion used to select these red sources. Dotted lines indicate 
the interstellar extinction path according to the extinction law 
of Mathis i ll 9901 1. The two nearly vertical dotted lines mark ex- 
tinctions of Ay = 10 and 20 mag. Note that the colour scales are 
in instrumental units. 
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Fig. 3. Colour-colour diagram of the IRAS 2 1 346+57 14 glob- 
ule. Symbols as in Fig. [2] 



In a first step, we executed several simulations with stan- 
dard Galactic optical extinction (0.7 mag/kpc) and found that 
the number of objects with H-K > 0.8 mag (and even with H- 
K > 0.5 mag) is zero. This Galaxy model, however, does not in- 
clude ultracool dwarfs with spectral types L and T. According 



to the DUSTY models of Chabrier et al. ( I2000l >. evolved ul- 
tracool objects with T c ff = 1700 K have an H-K colour above 
0.8 mag, and thus could be selected as red objects with our 
diagrams. On the other hand, objects with T c ff below 1300K 
are too faint to be seen in our images. For this range of ef- 



fective temperatures, Gizis et al. (I2000> give a space density 
of 2.11 • 10~ 3 pc~ 3 . Our limiting magnitude constrains the 
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Fig. 4. Colour-colour diagram of the IRAS 21352+5715 glob- 
ule. Symbols as in Fig.|2 




0.5 1 
H-K (mag) 



Fig. 5. Colour-colour diagram of the IRAS 21445+5712 glob- 
ule. Symbols as in Fig.|2] 



distance at which we could detect these ultracool objects to 
< 50 pc. With these values, the number of L-type objects in a 
1 sqare degree field is rs 0.03. From this estimate, we conclude 
that the number of Galactic L-type stars in our red sample is 
negligible. 

For this first calculation, however, we made the assump- 
tion of uniform Galactic extinction, which is not given in our 
fields, since we observe regions with high extinction in the di- 
rection of IC 1396. The Besangon Galaxy model is able to in- 
clude clouds with a given distance and extinction. In a second 
set of simulations, we therefore assumed a cloud at a distance 
of IC 1396 (750 pc) with various extinction values ranging from 
Ay = 3 mag to Av = 10 mag. The basic effect of such an addi- 
tional cloud is that the main sequence in the colour-colour dia- 
gram is split. While the foreground objects remain unreddened, 
the background stars are shifted along the reddening path. From 
these simulations, we found that the background stars begin to 
exceed the H-K = 0.8 mag limit if the extinction of the cloud is 
higher than Ay = 5 mag. Thus, we expect a significant number 
of reddened background stars for globules with Ay > 5 mag. 
According to the Galaxy model, these reddened objects will be 
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Fig. 6. Colour-colour diagram of the IRAS 2205 1+5848 glob- 
ule. Symbols as in Fig.|5] For this globule the selection criterion 
for the reddened objects was shifted to H-K > 0.7 mag, because 
all objects in this field show reddening of 0.2 mag. 

mostly M type giants for a cloud with Ay = 5 mag, but with 
stronger extinction stars with earlier spectral types also would 
appear to have H-K > 0.8 mag. 

Based on these results, we can make the following conclu- 
sions about the nature of the reddened objects: We can exclude 
that our reddened objects contain a significant fraction of fore- 
ground stars. This includes also ultracool dwarfs with spec- 
tral types L and T. For globules with high extinction, we ex- 
pect a significant number of highly reddened background stars, 
mostly late-type giants, along the reddening path. The remain- 
ing reddened sources, particularly the objects whose colours 
place them below the reddening path in the colour-colour di- 
agram, should be young stellar objects (YSO) associated with 
the globules in IC 1396, which are intrinsically red because of 
excess emission from circumstellar matter. From photometry 
alone, it is not possible to unambiguously distinguish between 
YSOs and background stars. 

Five of our nine globules clearly show a cumulation of red- 
dened objects, whereas the number of reddened objects in the 
remaining four globules is <10. It is now interesting to anal- 
yse whether an excessive number of red objects is caused by a 
high density of YSOs in the globule or a high cloud extinction, 
leading to a high number of background stars which appear red- 
dened. From the position of the reddened objects in the colour- 
colour diagram, we are able to get a rough idea of whether they 
are mostly YSOs or background stars, since background stars 
will be reddened along the reddening path, whereas YSOs may 
also appear below the extinction path. We define as YSO can- 
didates all objects whose photometry does not rule out the pos- 
sibility that they are located below the reddening band. This in- 
cludes all objects with photometry in all three bands, which are 
below the reddening band, but also objects with upper limit in 
J, that could be located below the extinction band, and all red- 
dened objects for which only K-band magnitudes are available. 
For each globule, we counted the number of YSO candidates 
and the total number of reddened objects. For globules with- 
out a large reddened population, it was found that only 77 % 



of the red sources are YSO candidates, whereas the remain- 
ing objects are in the extinction path. On the other hand, for 
globules with a large reddened population, the fraction of YSO 
candidates is 68%, i.e. significantly lower than in the 'empty' 
globules. The most straightforward interpretation of this result 
is that the globules with many red objects have both a high den- 
sity of YSO and strong extinction, leading to a high number of 
reddened background stars. The other globules possess smaller 
extinction values and fewer young stellar objects. 

Assuming a distance of 750 pc for IC 1396 and no signifi- 
cant extinction, our reddened targets have absolute J-band mag- 
nitudes between 3.1 and 7.3 mag. This constrains the masses 
for these objects to a range roughly between 0.05 and 0.9 M Q 
(Baraffe et al. IT9981 . Note that there are also no brighter 
sources significantly below the reddening path but with H- 
K < 0.5 mag. Hence no unextincted Herbig AeBe stars are 
present. Since the globules show J-band extinction of at least 
0.9 mag (see Tabled, there could be a few higher mass stars in 
our sample. On the other hand, it is unlikely that our faintest 
objects have substellar masses, rather than just being strongly 
influenced by extinction. 

4. Extinction maps 

4.1. Method 

We determined extinction maps of the region 2l'?3 to 22 l }2 in 
R.A. and 56° to 60° in DEC (J2000) using accumulated star 
counts (Wolf-diagrams) in the 2MASS database. An example 
Wolf-diagram is shown in Fig. IC.ll Stars are counted in 3' x 3' 
boxes every 20" down to the completness limit of the cata- 
logue. A co-centred 1 ° x 1 ° sized field was chosen as "unex- 
tincted" comparison. Minimum and maximum traceable ex- 
tinction values are calculated and shown in Fig. IC.2l depend- 
ing on the resolution used. We can trace about one to 20 mag 
optical extinction. Details of the star count method and the de- 
termination of the detection limits can be found in Annendixlcl 

4.2. Results 

The extinction map of the whole IC 1396 field shows several 
regions of high extinction (see Fig.[0. These are the globules 
around the star HD 206267, and a further group towards the 
north-east. We defined a 'globule' as an object with an ex- 
tinction of at least 3 a above the noise level and a size of at 
least 9 square arcminutes. With these criteria, an automated 
search using the SExtractor (Bertin & Arnouts 1996 1 revealed 
20 globules in IC 1396. Five of these globules were contained 
in the target list for our near-infrared survey with MAGIC from 
Schwartz et al. J199 li . The globules in this list that could not 
be detected in our extinction maps possess sizes and masses 
below our detection limits (see also Patel et al. 1995 >. Adding 
the remaining five objects from this list, there are 25 dark 
clouds in the field. The large region of high extinction at about 
(2l'. 1 73:58.8°) is caused by the bright star /iCep, which pre- 
vented the detection of stars in its vicinity. Similarly the star 
CCep at (22^19:58.2°) causes such a fake globule. The full list 
of all globules with their coordinates and known identifications 
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is given in Tabled Figure^ shows as examples the extinction 
maps for the globules where we detected a large number of 
reddened sources. 

To measure the mass we integrated the total extinction in 
the globules (Ay*). Only regions with an extinction above the 
one sigma noise level where considered, and the outer radius 
of the globule was taken as the size given by the SExtractor 
software (FWHM, assuming a Gaussian core). From these 
integrated extinction values we determined the mass of the 
globules using the fact that the column density N(H) of hy- 
drogen atoms can be expressed as 6.83 10 21 cm~ 2 x Ay/Ry 
where Ry is typically 3.0 (see e.g. Mafhis 119901 . Using a 
distance of 750 pc for the globules, and the 20" pixel scale 
in our maps, the total mass in a globule can be determined 
from the integrated optical extinction (Ay*) in the globule by: 
Mgiob [M ] = 0.098 Ay * [mag]. Note that the obtained extinc- 
tion values are averaged over the box-size of 3' x 3'. If the dust 
is concentrated in a smaller area, the obtained masses are lower 
limits. The error of the estimated mass depends on the mean 
extinction within the globule compared to the noise level in 
our map. As a typical value we find that the mean extinction is 
about 5 u of the noise level and hence the uncertainties in the 
inferred masses are about 20%. In the case of IC 1396N we 
can compare our mass estimate with literature values. Wilking 
et al. J1993> estimate 380 M Q and Serabyn et al. dl993> give 
480±120M Q within the central 0.3 pc. Given the large un- 
certainties and the fact that the object is rather small and the 
mass is concentrated close to the actual star forming core (see 
e.g. Codella et al. 2001), our lower limit for the whole glob- 
ule of 300 M Q agrees well with these estimates. Weikard et 
al. ( 1996 > give an estimate of 12000M© for the total mass 
of their mapped region (about six square degrees centred on 
HD 206267). The mass estimated from our extinction maps in 
the same field is 9000 M©, reasonably close. 

Within the errors of about 20 % the total optical extinc- 
tion values and globule masses, obtained from the J- and H- 
band 2MASS data, are consistent. Some of the masses calcu- 
lated with the K-band data differ by a larger amount. There 
are two groups of objects: 1) The mass from the K-band data 
is much smaller than the mass obtained from the J- and H- 
band data. The reason could be that the globule contains an 
embedded cluster of stars. These stars are easier to detect in K 
and hence decrease the apparent extinction. In these cases the 
given masses Mj.h are lower limits for the globule mass. 2) 
The mass from the K-band data is much larger than the mass 
obtained from the J- and H-band data. An explanation could be 
that the conversion of extinction in the K-band to optical extinc- 
tion cannot be performed following Mathis et al. (1990) that 
assume an opacity index of 1.7. A lower value might be valid 
in these globules. Note that this effect could counter-balance 1) 
in the case of an embedded cluster. 

5. Globule properties 

With the results of Sect. [5] and |4] a homogeneous sample of 
globules is now at our disposal (see Table [0- To evaluate pos- 
sible trends in this dataset, we selected the following glob- 
ule properties and tested them for correlations: projected dis- 



tance of the globule from the exciting O star HD 206267, pro- 
jected size of the globule, number of YSO candidates within 
the field of the globule, peak extinction value, total mass, mean 
column density within the globule (total mass divided by the 
size). For each pair of these properties, we fitted their relation- 
ship linearly and calculated the correlation coefficient r. Since 

the value t = r^J ^S r V - (N: number of datapoints) follows 
Student's i-distribution, this gives us the probability that the 
two parameters show indeed a significant correlation. In or- 
der to perform these tests, we set globule sizes to 30 square 
arcminutes and masses to 120 M©, if they could not be deter- 
mined. These values are the upper limits from our detection 
procedure for globules in the extinction maps. As peak extinc- 
tion, we used the values from the J-band, since they have the 
best signal-to-noise ratio. 

The best correlations are obtained when comparing size and 
mass, peak extinction and mass, and size and peak extinction. 
The false alarm probabilities (FAP) for linear relationships be- 
tween these three parameters are below 10~ 3 %. These are ex- 
pected correlations, since the mass is proportional to the mean 
extinction and the size of the globule. The correlation of size 
and peak extinction might be due to our method that smooths 
the extinction over a 3' x 3' field, naturally leading to lower 
peak extinction values for small globules. 

Additionally, we found a non-obvious correlation between 
mass M and distance d from the exciting star: Globules far- 
ther away from HD 206267 have on average larger masses. The 
least-square fit gives M[M e ] = 8.6 d [pc] +84, and the FAP for 
this correlation is 0.078 %. We repeated this test excluding the 
four globules with d > 40 pc, because it might be that these re- 
gions are too far away from HD 206267 to be influenced by its 
radiation. Without the data points for these globules, the FAP is 
still only 2.2 %. Thus, there is a significant connection between 
globule mass and distance from the exciting star (see Fig. [8}. 

There are two possible explanations for this correlation: 1) 
The mass loss rate of a cloud induced by photo-ionisation of 
the external layers is inversely proportional to the distance from 
the exciting star (see e.g Codella et al. 1200 It . Hence, more dis- 
tant globules do not lose as much mass due to evaporation. 2) 
Globules closer to the exciting star contain more young objects 
(see below) and hence the mass estimates of these globules lead 
to too low values due to the star count method. Connected to the 
distance-mass correlation are correlations of distance and glob- 
ule size as well as distance and peak extinction value. These are 
consequences of the interrelations described in the above para- 
graph. 

The remaining combinations of parameters show no signif- 
icant correlations, i.e. with FAP below 10%. Particularly, we 
see no correlation between distance and number of YSO can- 
didates. This can be used to rule out the second explanation 
for the mass-distance correlation. This relationship is therefore 
most likely related to mass loss via photo-evaporation, as ex- 
plained above. For all correlation tests that use the number of 
YSO candidates, we are restricted to those globules in which 
red objects were identified (see TableQ. Thus, small number 
statistics hamper the correlation tests for these parameters. 
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Fig. 7. Contour plots of the extinction obtained from the J-band 2MASS data for the globules 1 1 (upper left), 3 (upper right), 
9 (lower left), and 4 and 5 (lower right). Contours start at 1 mag optical extinction and increase by 1 mag. Conversion from 



extinction in the J-band to optical extinction was done according to Mathis et al. ( 1990). Dashed lines mark the borders of the 
observed NIR fields. + signs indicate reddened objects in the reddening path, while circles mark the YSO candidates, and arrows 
the principle direction towards the 06.5 V star HD 206267. Note that the contours are slightly smoothed and hence peak extinction 
values are somewhat lower than in Tabled 



How do these findings fit in a scenario where the radiation 
pressure of the O star is the primary triggering mechanism for 
star formation in the globules? In a first order interpretation, we 
assume that the density of young objects S/A 3 / 2 (S - numer of 
YSO candidates; A - projected size of the globule) is positively 
correlated with the mass of the globule M and the pressure P 
exerted by the radiation flux $ from the star: S/A 3 / 2 oc M x P. 
The radiation pressure or ionising flux is inversely proportional 
to the square of the distance (P oc $ oc d~ 2 ; Codella et al. 
I200T1 and the mass is proportional to the distance (M oc d as 
shown in Fig. [8}. Hence we should expect a correlation of the 
star density with the inverse distance from the star: S/A 3 ' 2 oc 
M x P oc Indeed, such a correlation is seen tentatively 
in our data. A linear fit gives S/A 3 ' 2 = Q.555d-\ with a FAP 
of 2.5 % (see Fig.[9j. If we exclude the one globule with d > 
40 pc, the FAP increases to 5.0%. Although this result clearly 
needs to be substantiated with more datapoints, it tentatively 
confirms our initial assumption that the star forming activity is 
driven by the radiation pressure of the O star. 

We re-examined the conclusion of Sect. 13.21 where we 
found that globules with a rich population of young stars 
also show high extinction leading to a high number of red- 



dened background objects. Although there is no correlation be- 
tween extinction Aj and the number of YSO candidates N, 
the data show a clear tendency: All globules with N < 10 have 
Aj < 1.1 mag. On the other hand, 50% of the globules with 
N > 10 exhibit Aj > 2 mag. Thus, globules that harbour many 
red objects tend to show high extinction values, confirming our 
result from Sect. 13.21 



We further investigated the positions of the YSO candi- 
dates in the globules with respect to the direction of the 06.5 V 
star HD 206267. Figure^] shows the positions of the reddened 
sources overplotted on the extinction map for the five globules 
where a large number of reddened objects was found (the two 
overlapping fields are combined in one figure). In three of the 
fields the YSO candidates (circles) seem to be preferentially 
positioned towards the direction of the exciting 06.5V star. 
Objects within the reddening path (plus signs) seem to be more 
concentrated towards the high extinction regions. This picture 
agrees well with the expectations of triggered star formation 
via radiation-driven implosion, although this behaviour is not 
present in all investigated globules. 
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Table 1. All detected globules in the IC 1396 region. The upper part of the table lists the objects observed with MAGIC in IHK, 
while in the lower part the additional globules detected in our extinction maps are given. Column 4 lists the projected distance 
of the globule from the exciting star HD 206267. The size in square arcminutes is the value we obtained from the SExtractor 
software (a ? marks the globules that were not detected in the extinction maps.). Column 6 contains the number of reddened 
objects (in brackets: number of YSO candidates as defined in Sect. l3.2t . In columns 7 to 9 we give the peak extinction values in 
the three filters obtained from our extinction maps, while columns 10 to 12 list the masses estimated for the globules. * The given 
masses are measured for both globules together. 
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6. New HH objects and H 2 outflows 

In our observed field there are several known Herbig-Haro ob- 
jects and outflows. These are mainly the well investigated flows 
in the IC 1396N globule (HH 589-595 found by Ogura et al. 
l2002l and HH 777-780 discovered by Reipurth et al. l2003l The 
H 2 outflow from IC 1396 W which was shown in Froebrich & 
Scholz ( 2003) has no published optical emission counterpart. 
Connected to IRAS 21388+5622 is the HH588 object (Ogura 
et al. l2002l l. In the IRAS 22051+5848 field a giant outflow 
HH 354 was found by Reipurth & Bally J200T1 >. 

In our [SII] images we were able to re-discover all of the 
known HH objects around IC 1396N. The whole field is full 
of extended emission line objects (visible as excess emission 
in the [SII] images compared to the I-band fluxes). These fil- 
aments are excited by the strong UV emission of the 06.5V 
star HD 206267. Many of these filaments show a similar ap- 
pearance to some of the HH objects near IC 1396N. Hence, it 
is very difficult to decide from the shape of the emission alone 



Table 2. Coordinates of the brightest knots in the newly de- 
tected Herbig-Haro objects. * Counterparts of the optical knots 
detected in H2 in Froebrich & Scholz J2003> . 

Object H2 q(I2000) 5(J2000) 
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if we see an outflow or just UV excited filaments. This might 
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Fig. 8. Globule masses vs. distance from the exciting O star. 
A linear least-square fit is shown as dashed line. This corre- 
lation has a FAP of 0.078 %. If we exclude the four globules 
with d > 40 pc, the correlation is still significant with a FAP of 
2.2 %. 
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Fig. 9. Diagram of star density vs. inverse distance from the 
exciting O star. A linear least square fit is shown as dashed line. 
The correlation has a FAP of 2.5 %. If we exclude globules with 
d > 40 pc a FAP of 5.0 % is still reached. 



be the reason why HH777 was not mentioned by Ogura et al. 



(2002 ) although it is clearly visible in their image. 

We also detect the HH354 and HH588 flows from 
IRAS 22051+5848 and IRAS 21388+5622, respectively. There 
is an additional emission feature 1.'4 east of HH 588 NE2, con- 
sisting of two fuzzy blobs (for positions see Table[2j. This ob- 



ject was slightly outside the field shown in Ogura et al. y002 1. 
The knot-like appearance of this object suggests an HH object 
(HH 588 NE3) connected to the HH 588 flow, rather than being 
a UV-excited emission feature (see Fig.llOi. None of the optical 
emission features of HH 588 is detected in our H2 images. 

We detected two new groups of HH objects in our op- 
tical field. We found the H 2 flow from IC1396W also in 
[SII] and call it HH 864 (s ee Fig-EJ- The bright H 2 knot 2-j 
(Froebrich & Scholz l2003l has a [SII] counterpart (HH 864 A) 
which shows a bow shape with two maxima in the emission. 
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Fig. 10. [SII] image of the HH588 region. The newly detected 
emission knots HH 588 NE3 are seen in the upper left corner. 
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Fig. 11. Our [SII] image of the outflow from IC1396W 
(HH 864) shows clearly some of the brightest emission knots, 



previously detected in H 2 by Froebrich & Scholz J2003> . 

The two knots (2-b and 2-d) in the north-eastern lobe are vis- 
ible in the optical emission line also (HH864C). Knot 5 (50" 
east of 2-1) has an optical counterpart as well (HH864B). 
It is still not clear if this emission knot is connected to the 
IRAS 21246+5743 source or another (unknown) source in the 
IC 1 396 W globule. 

We find in our optical field a giant flow emerging from 
the IRAS 21445+5712 source, which is also called IC 1396E. 
There are two bright bow shocks heading to the north-west of 
the source, HH 865 A and HH 865 B (see Fig.[l2). The distance 
of about 0.2° of the terminating bow HH 865 A from the source 
makes this flow 2.6 pc in length in one lobe (assuming a dis- 
tance of 725 pc). No counterflow is found in our images. The 
globule was also observed in H 2 . The closer of the two bow 
shocks is unfortunately just outside our NIR field and hence 
could not be detected. There is no H 2 emission in our image 
that might be connected to this flow. Only a few very faint fea- 
tures are found south and east of the IRAS source. These are 
UV excited emission coinciding with cloud borders. 

In IRAS 21312+5736 we see some faint H 2 filaments 
that might be UV excited. The field of IRAS 21324+5716 
is full of bright H 2 emission features coinciding with the 
optical emission. This applies also for the neighbouring 
fields IRAS 21346+5714 and IRAS 21352+5715. No H 2 fea- 
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Fig. 12. The new flow HH 865 in the [SII] filter, heading away 
from the possible driving source IRAS 21445+5712. 



flow from IC 1396 E. Further a new emission knot belonging to 
the known HH 588 object was discovered. 
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tures are detected near IRAS 21354+5823, IRAS 21428+5802, 
IRAS 21539+5821 and IRAS 22051+5848. 



7. Summary and discussion 

NIR observations of IC 1396 in conjunction with extinction 
maps obtained from 2MASS data reveal star forming activ- 
ity and a large number of globules in this region. Twenty five 
globules were identified using our extinction maps and the list 
of Schwartz et al. (1199 II . Four of them were previously un- 
catalogued in the SIMBAD database. In all but four cases the 
masses (or at least lower limits) of the globules could be de- 
termined. Also the size could be measured properly for all but 
seven objects. 

For nine globules observed deeper than 2MASS in this 
work, the content of heavily reddened objects was derived by 
means of colour-colour diagrams. Five of these globules exhibit 
a rich population of red objects. At least half of these objects 
are good candidates for young stellar objects, the remaining 
half is probably contaminated by reddened background stars. 
The five globules with many red objects include the targets 
with the highest extinction values, suggesting a correlation of 
the strength of the star formation activity with the mass of the 
globule. 

Star formation in small globules is often thought to be 
strongly influenced by the radiation pressure of a nearby bright 
star. It was therefore investigated how the globule properties in 
IC1396 depend on the distance from the O star HD 206267. The 
masses of the globules show a clear positive correlation with 
the distance from this star. We conclude that evaporation due 
to photo-ionisation affects the mass distribution of the globules 
around HD 206267. Our data are consistent with a scenario in 
which the radiation pressure from the O star regulates the star 
forming activity (expressed as density of young sources) in the 
globules, in the sense that the radiation pressure compresses the 
gas and thus leads to enhanced star formation. 

Our optical data lead to the discovery of several new HH 
objects. These are the counterpart of the previously known H2 
emission of the flow from IC 1396W and a new parsec scale 



Appendix A: Reliability of instrumental 
photometry 

Several tests were executed to evaluate the reliability of the 
instrumental photometry. Since we observed each globule at 
least twice (and six of them three times), we can compare the 
instrumental magnitudes from different mosaics of the same 
region to look for systematic mismatches. With one excep- 
tion (see below), we found good agreement within ±0.05 mag. 
The photometry from mosaics which were taken under non- 
photometric conditions (see Sect. 12. li shows the highest devia- 
tions, as expected. These deviations, however, affect all broad- 
band mosaics equally, thus they are not visible in the colours. 
For one field (IRAS 21539+5821) we see large systematic pho- 
tometry mismatches between all mosaics, which are also vis- 
ible in the colours. The most probable explanation is variable 
weather conditions, e.g. cirrus clouds which we did not recog- 
nise during the observations. We excluded this field from all 
further analysis. 

As mentioned in Sect. 12.11 all broad-band images were ob- 
tained at low airmass; the maximum airmass difference is 0.45. 
Since the NIR extinction coefficient on Calar Alto is usually 
below 0. 1 mag/airmass (Hopp & Fernandez l2002t . the differen- 
tial extinction offsets between the instrumental magnitudes of 
different fields are safely below 0.05 mag. More important, the 
extinction coefficient in the NIR is nearly wavelength indepen- 
dent. Therefore, differential extinction does not significantly 
affect the colours of our targets. Thus, we did not perform an 
extinction correction. Another test of our photometry was en- 
abled by field overlap of the mosaics around IRAS 21346+5714 
and IRAS 21352+5715. Again we found no significant offsets 
between the colours of the targets which were detected in both 
fields. For these reasons, we conclude that the instrumental 
colours of the targets observed under photometric conditions 
are reliable. 

If available, the deep, stacked mosaics of the globules were 
calibrated by measuring the magnitude offsets to a mosaic ob- 
tained under photometric conditions and applying these offsets 
to the instrumental magnitudes of the deep mosaics. Thus, the 
photometry of the deep mosaics is now directly comparable 
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with all other data obtained under photometric conditions, in 
particular with the standard star main sequence. 

Appendix B: Colour-colour diagrams and 
absolute calibration 

For the colour-colour diagrams, we used the photometry from 
the deep mosaics if more than one mosaic is available, other- 
wise the photometry from the single mosaic. We only consider 
objects with errors below 0.2 mag. There are, however, only 
very few objects (typically below 1%) with larger errors. 

The bulk of the IC 1396 objects in the diagrams coincides 
with the late-type end of the main sequence. There is one ex- 
ception - Fig.[6]- which will be discussed below. This is ex- 
pected because most stars are late-type. Moreover, giant stars 
also concentrate towards the late-type end of the main se- 
quence. Since our main sequence ends at spectral type K8, it 
is not surprising that there are many field objects above the up- 
per end of this sequence. These could be M-type dwarf or giant 
stars. The bulk of the field objects is (again with one exception) 
always in the same position in each diagram, which makes us 
confident that the colours from the different fields can be com- 
pared with each other. 

There is one field (IRAS 22051+5848, see Fig.®, where 
the bulk of the datapoints is clearly offset by about 0.2 mag 
in (H-K) and 0.2 mag in (J-H), i.e. roughly in the direction of 
the reddening vector. This field shows large-scale structures 
of nebulosity and large voids without any objects, indicative 
of strong extinction. Thus, in this field background objects 
are probably significantly reddened because their light must 
pass through extended regions of dust. The 0.2 mag offset in 
both colours then corresponds to a mean visual extinction of 
A v ~ 2 mag. 

We are now interested in detecting objects whose position 
in the colour-colour diagram indicates significant intrinsic red- 
dening, i.e whose position is clearly shifted in the direction of 
the extinction path or even below these lines. The following 
criteria were used to select such reddened objects: 

a) If an object is detected in all three filters, its H- 
K colour should exceed 0.5 mag: H-K > 0.5 mag. (For the 
IRAS 22051+5848 field, we require H-K > 0.7 mag, because 
all objects in this field show a 0.2 mag shift in H-K, see above.) 
Objects that satisfy this condition are shown as filled squares 
in the diagrams. 

b) If the object is only detected in H and K, 
we demand again that H-K > 0.5 mag (or > 0.7 mag for 
IRAS 22051+5848). These objects are marked with an 'arrow 
up', and their J-H colour is a lower limit estimated from the 
H-band photometry and our sensitivity limit in the J-band. 

c) If the object is only detected in K, we can only deter- 
mine a lower limit for the H-K colour by subtracting the K- 
band photometry from the sensitivity limit in the H-band. We 
require that this lower limit is > 0.5 mag (or > 0.7 mag for 
IRAS 22051+5848). These objects are shown with an 'arrow 
to the right', and their J-H colours are chosen arbitrarily. 

We selected all objects that satisfy one of these condi- 
tions and examined these targets in the images. We rejected 
all objects that are clearly not star-shaped (i.e. galaxies or H 2 



emission knots), as well as spurious detections (i.e. cosmics or 
spikes of bright stars). The remaining objects are good candi- 
dates for objects with intrinsic reddening in the globules. This 
selection might be incomplete, since some young objects could 
have H-K < 0.5, and thus would fall out of our selection crite- 
ria. As noted above, the colour-colour diagrams contain only 
objects with errors below 0.2 mag. With larger error bars a re- 
liable separation between reddened and unreddened objects is 
not possible. In most cases, however, the objects with errors 
> 0.2 mag do not appear strongly reddened. 

After this selection process, the globules clearly fall in 
two groups: Four of them show only very few reddened ob- 
jects, i.e. their number is < 10. The remaining five glob- 
ules harbour more than 20 reddened objects, in two cases 
(IRAS 21346+5714, IRAS 22051+5848) the number of red- 
dened objects is larger than 50. For these five globules we show 
the colour-colour diagrams in Figs. 12161 

As noted above, the colours in the diagrams are instrumen- 
tal. We determined, however, the offsets to the photometric sys- 
tem of the 2MASS catalogue, which is available online 3 . For all 
detected objects, we searched for counterparts in the 2MASS 
database. In each field, several hundred common objects were 
identified. For these objects, we determined the average offset 
A between 2MASS photometry and our instrumental magni- 
tudes in J, H, and K. These offsets (called A,/, Ah, Ak in the 
following) include a zero-point as well as extinction correction. 

The offsets in a certain band are not constant for all glob- 
ules; they differ by as much as 0.3 mag, depending on the 
weather conditions during the observations (see above). The 
difference between the offsets, in particular the values of A,/ — 
Ah and Ah — Ak, are comparable for all globules, con- 
firming again that the relative colours are reliable. We obtain 
Aj - A H = 0.03 ± 0.06 and A H - A K = 0.29 ± 0.08. Thus, 
our colour-colour diagrams can be roughly transformed into 
the 2MASS photometric system by adding 0.3 mag to the H- 
K values. We note, however, that this transformation does not 
include colour terms, which means that the offsets are proba- 
bly not constant for all colours, in particular not for the main 
sequence and the reddened targets (see Froebrich & Scholz 
2003)). This was the main reason why we kept the diagrams 
in the instrumental colours, as described above. 

The absolute calibration also allows us to estimate the sen- 
sitivity limit of our survey, i.e. the faintest objects for which 
a 5a detection is possible. We reach about 17.0 mag in J, 
16.5 mag in H, and 16.0 mag in K. For comparison, the 2MASS 
catalogue for this region contains objects down to 16.9 mag in 
J, 16.1 mag in H, and 15.6 mag in K. Our images are thus sig- 
nificantly deeper than the 2MASS catalogue for this region, 
particularly in the H- and K-band. 

Appendix C: Extinction map determination 

A determination of absolute magnitudes of the stars within the 
globules, as well as an estimation of the dust mass in the glob- 
ules, requires a measurement of the extinction. As discussed 
above, the assumption of a uniform Galactic extinction is not 

3 see http://www.ipac.caltech.edu/2mass 
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Fig. C.l. Example of a Wolf diagram for the three filters (J - 
solid; H - dashed; K - dotted) at a = 21 h 52 m 00 s 5 = 57°00'00" 
(J2000). The accumulated number of stars N is averaged over 
a 1° x 1 "field and normalised to an area of nine square arcmin- 
utes. The vertical lines correspond to the completeness limit in 
JHK (line type as for the Wolf diagrams) at which we measured 
the extinction. The mean slope <X> varies with the filter (0.36 
for J and 0.34 for H and K) and also with the position on the 
sky. As can be seen the slope is also not constant, due to ex- 
tinction from the ISM (the slope in J ranges from 0.38 (at 1 1 to 
12 mag) to 0.30 (at 15 to 16 mag)). 

valid in our field. There are local extinction enhancements con- 
nected to the globules and even within these small clouds the 
dust and so the extinction is not uniformly distributed. Hence 
we need to determine the local extinction enhancements for 
our globules, compared to the neighbouring star field. These 
extinction enhancements will be used to estimate the mass of 
the globules by adopting a distance of 750 pc. Measuring the 
extinction from colour-colour diagrams is relatively difficult 
due to the photometric errors and the large scatter of the main 
sequence (see e.g. our Figs. 12161 . Hence this method only al- 
lows a rough estimate of at best ± 3 mag of the optical ex- 
tinction averaged over the whole globule and is in particular 
not able to determine the extinction/mass distribution within 
the globules. A more accurate method is to create accumula- 
tive star counts (Wolf diagrams) as described e.g. in Kiss et 



al. (|2p00}. Since we expect high extinction values within the 
globules (Ay > 5 mag), creating these diagrams from NIR ob- 
servations obviously is the best choise. The 2MASS database 
provides an ideal basis for this purpose. 

Using star counts to estimate the extinction is based on sev- 
eral assumptions: 1) The stars are equally distributed and all 
apparent voids or less densely populated regions are caused 
by extinction. 2) All stars possess the same absolute bright- 
ness. 3) The completeness limit of the input catalogue does not 
depend on the position in the sky. Since in general these as- 
sumptions are not valid, we performed our star counts in a way 
that ensures as small as possible errors due to the assumptions 
made. 1) As reference or control field we choose a 1° x 1° field 
around the position where the stars are counted. This running 
average ensures that large scale variations in the average star 
numbers due to Galactic position or large scale structures are 




Resolution ['] 



Fig. C.2. Minimum and maximum values of optical extinction 
that could be traced with our extinction maps obtained from the 
2MASS data in JHK, using different box sizes. Ay values are 
determined from extinction in JHK using the conversion factors 



given in Mathis et al. ([T990i. 

corrected. 2) The box size in which we count the stars for each 
position was varied between l'x 1' and 5'x5'. In boxes larger 
than 3'x3' of the order of 100 stars are enclosed on average. 
This ensures that we enclose a 'typical' sample of stars in our 
box, that represents the 'average' star very well. 3) The com- 
pleteness limit of the 2MASS catalogue varies over the field. 
Star counts where hence only performed down to a magnitude 
to which the catalogue is complete over the whole field (J=16.6, 
H=15.6,K=15.2). 

In order to investigate not only the extinction distribution in 
our globules but also the large scale dust distribution in the en- 
tire Cep OB 2 region, we created an extiction map of the whole 
field, ranging from 21^3 to 22^2 in RA. and from 56° to 60° 
in DEC (J2000). In Fig. ^ we show an example of the extinc- 
tion map obtained from the J-band 2MASS data. Since we had 
a large amount of data to process we decided to parallelise the 
problem. We wrote an MPI-parallel C++ code in order to get 
maximum performance from our 32-node hyperthreaded clus- 
ter of the so-called "Beowulf type. It took 13 hours to process 
all the data on this cluster - without the parallelisation of the 
problem it would have taken of the order of 60 times as long 
(about 33 days). 

For the star counts we selected all 2MASS objects in our 
field with a signal-to-noise ratio of larger than five (quality flag 
A, B, or C). With this selection criteria we extracted 841908, 
783474, and 635478 objects in JHK, respectively. This con- 
verts to an average density of stars of about 1 star per 20" x 20" 
field. Hence, we performed the star counts every 20" in or- 
der to obtain the final extinction map. This method ensures 
a maximum gain of information from the 2MASS catalogue. 
Figure lcTl shows an example of a Wolf diagram in each of 
the three filters at a = 21 h 52 m 00 s <5 = 57°00'00" (J2000). In 
order to find the best compromise between box-size and ex- 
tinction regime that can be traced, the box-sizes were varied in 
steps of 0.'25. Depending on the box-size the average number 
of stars in the boxes varies and hence we can trace different 
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extinction regimes. In Fig. IC.2l we plot the minimum and max- 
imum value of traceable extinction (converted to Ay) with our 
method in IC 1396, depending on the chosen box-size and the 
three 2MASS filters. The lower limit (A™ in ) is calculated from 
the three sigma noise level of the determined extinction in the 
whole Cep OB 2 field. This limit can be transformed into a ratio 
F^/Fback by 

Af n = ~*logfl-=^_). (C.l) 

A \ r back / 

F a is the noise of the number of stars in the star count map. If 
this number of stars is missing in the presence of Fb ac k stars 
in the comparison field, an apparent extinction of A™ ln is de- 
tected. The factor X is the slope in the Wolf-diagram. Note that 
this factor, as well as the number of background stars varies 
with the position and filter. The maximum tracable extinction 
A™ can now be determined by assuming that only F CT stars 
are present. 

= --* log ( J (C.2) 

X \ Fback / 

This is a conservative assumption, because F^ represents ac- 
tually the noise of the number of background and foreground 
stars, but only the latter determines the maximum of tracable 
extinction. Note: If we assume that the scatter of foreground 
stars accounts only for half of F CT , the maximum tracable ex- 
tinction values increase by about 4, 6, and 10 mag optical ex- 
tinction for JHK, respectively (assuming X = 0.3). We selected 
a box-size of 3' x 3' to perform all our measurements. The 
lower limit (A" lm ) is also comparable to the error of our de- 
termined extinction values. As can be seen in Fig. IC. 21 this is 
much lower than the estimated accuracy of ± 3 mag when the 
extinction is inferred from the colour-colour diagrams. 
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